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FATIGUE ANALYSIS OF COMPOSITE MATERIALS USING THE FAIL-SAFE CONCEPT 


by G. Stievenard 

Aerospatiale Helicopter Division 
0. INTRODUCTION 

The helicopter by design presents numerous sources of vibratory 
excitation and as a result many components are under stress loads. 
The main concern of the manufacturer is to produce components that 
meet highly strict safety requirements while remaining competitive 
in terms of mass and cost-price considerations. A fatigue regula- 
tion established jointly between SNIAS and Official French Services 
requires that the individual failure probability of a helicopter 
component must not exceed 10 ^ during a given lifespan. 

The development of composite materials over the pest fifteen 
years has made it possible for Aerospatiale to design new rotors 
that meet the specifications of the abovementioned regulation. 

These materials have a good fatigue life with evident fail- 
safe cnaracteris tics (easily detectable onset of delamination and 
slow propagation speed) . 

TJiese properties have made it possible to develop metliods of 
analysis without increasing tlie risk of failure during flight. 

If we let R1 be the probability of having a craclc on a flight 
component and R2 be the probability of seeing this crack propagate 
between two scheduled inspections, the global failure risk is the 
product of R1 X R2 and it is this product that must not exceed iO 

The objective of tliis report is to describe the method used by 
SNIAS permitting the calculation of these two risks R1 and R2. 
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This method is based on a full study of fatigue characteris- 
tics and crack propagation rates in the material used. 

1. CALCULATION OF RISK R1 (Risk of Having a Crack During Operation) 

The first stage of the study consists of determining the equa- 
tion of the S/N curve of the material used to manufacture the heli- 
copter component. These equations are generally expressed as shown 
in figure 1 below: 



From the tests of samples performed in fairly large numbers, 

:t is possible to calculate the coefficients A and alplia which are 
characteristic of a material. The metiiod used is a linear regression 
along log N as a function of set S (fig. 2) ; 



Fig. 2 
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The second calculation hypothesis consists of assuming that 
there are S/N curves correlated to the mean curve of the same 
equation A which follow a normal law of logarithmic probability. 

These curves are therefore curves of crack isoprobability 
(fig. 3): 



Fig, 3 

This hypothesis which may seem arbitrary at first is well 
verified by tlie experiment and statistical testing. 

The experiment that we have just described would be too costly 
if it were performed entirely on aircraft components. Also, the equa- 
tion of the S/N curve is determined on specimens and is considered 
to be applicable to only a few specimens of the components we are 
testing (G mini pts) fig. 4: 
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Formulation of the Calculations 
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Key: a-Number of test points; b-risk under consideration; 

c Mean fatigue limit at 10^ cycles; d-standard devia 
tion; e-Normal law coefficient (5), 


kr ■ margin corresponding to risk r in the simple normal relationship, 
kq ■ margin corresponding to risk q in the simple normal relationship, 
q = coefficient of confidence. 

Knowing the 6 test results, it is now possible to calculate Am 
and the Ar's corresponding to different risks r. 

Operational characteristics of a helicopter are known, either 
by recording flight parameters during several hundres of hours or 
by information obtained from users. We thus have a spectrum giving 
for each flight configuration the percentage of time corresponding 
to it. For each of these configurations, the flight stresses are 
recorded using extensometric gages. 


These stresses associated with their percentage serve to 
calculate the lifespan by applying the MINI’.R relationship (see 

Table 1 



Key: a-% use; b-Flight stress; c-Fatigue limit at risk r; d-Number 

of allowable cycles; e-Number of flight hours N x t per Me 
(t : time for 1 Megacycle); f-Hovering; g-Vi transition; 
h-Level flight; i-Approach; j -Landing. Mean time between Tallurc 

Fig. 6 
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Fig. 7 

2 - CALCULATION OF RISK R2 (Risk of having a catastrophic failure 
between two scheduled inspections) 

The tests performed on components to calculate risk R1 are 
suspended at the first sign of damage and are resumed fnr fligiit 
loads using a program block and taking into account: 

-the percentage each configuration is used, 

-the chronological order of these configurations (in order to 
reproduce possible delay effects of the propagation due to changes 
in stress levels). 

(See load spectrum on Table II fig. 8). 

These tests were carried out to the failure point with measure- 
ments and recording of crack lengtns and resistance of the components. 

If we call t^£ the time elapse between time t. corresponding to 
the appearance of the first small damage and time t^ corresponding 
to catastrophic failure, it is possible to match a probability 
relationship to this random variable t^^ found during the test. 
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Table II 

Spectrum of Programmed Loads 
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2 31.5 
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1 
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81 
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83.5 

2 28 

100 
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100 
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224 
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64 

s 

84 

2 168 

100 
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176 
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61 
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320 
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Hovering; b- 

Take-off; c 

-Minimum 

cruise 

speed ; 


d-Turning; e-Lconomy cruise speed; f-Turning; g-Economy 
cruise speed; h-Turning; i-Maximum cruise speed; j -Turn- 
ing; k-Maximum cruise speed; 1-VNE; m-VNE turning; 
n-Maximum cruise speed; o-Hconomy cruise speed; p-Turn- 
ing; q-Minimum cruise speed; r-Turning; s-Speed reduc- 
tion; t- Autorotation ; u-Hovering; v- Approach; w-Rotor 
stop . 


Essentially, it seems that the normal logarithmic relationship 
fits the phenomenon the best. 


Since we know the law of probability, the mean value 
and the standard deviation, it is possible to calculate the 
allowable flight time as a function of risk (see fig. 9) 

The curve thus plotted corresponds to the risk of failure by 
considering that the first damage was noticed right after take-off. 
The curve of risk R2 should therefore be evaluated by assuming that 
this first damage may occur at any moment of flight. 
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Fig. 9 

If we let X be the flight hours between two inspections, we 
have two possible cases: 

(a) No damage appears between two given inspections (this 
case is not of interest to us since it carries over to case b) 
during the next inspections: 

(b) Damage does appear between two inspections. Only this 
case is of interest to us. We will therefore let the probability 
of damage appearing be equal to 1. I.e. tj^ is the instant in which 
the first noticeable damage appears (t counting from zero of the 
last inspection) and t 2 is the propagation time bringing about 
catastrophic damage to the component. It should be pointed out that 
tj and t^ are two independent variables. 

Probability j * < *-i < x ♦ <ix | ^ 

according to the preceding curve 

* cl 

Probability < x - x| • r ( x - x) 

where dP • ^ . r (x - x) 

X 

P(x) ■ 1 . / R (X - x) . dx ■ Rj 

X o' 

This is the mean value from 0 to X of the function 
explained previously (see fig. 10). 
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/ Failure probability by 
"1 assuming that the 1 st 
‘ damage occurs at the 
beginning of the flight. . 5 

X ^ 

1 

Failure probability by 
assuming that the 1 st 



damage occurs at any time 
during the flight. * c 


3 - CALCULATION OF RISK R (Or Risk of Failure During Flight) 

As we have alioady mentioned, this risk is found by multiplying 
R1 X R2 , the calculation method of whicii we have just described. 

R1 and R2 are presented in the form of two curves. The manu- 
facturer therefore has several possible choices, either to assign 
a long lifespan with frequent scliedulcd inspections of the component, 
or, conversely, a short lifespan with very spaced inspections. 

V/)iat is even better, these curves make it possible to modulate 
the inspection periods as a function of the component's flight 
time (f ig . 11 ) . 
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^Periods between inspections 
Fig. 11 



